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■ ABSTRACT 

^ , 

We have compiled a large sample of low-redshift active galactic nuclei (AGN) iden- 
tified via their emission line characteristics from the spectroscopic data of the Sloan 
^ . Digital Sky Survey. Since emission lines are often contaminated by stellar absorption 

On , lines, we developed an objective and efficient method of subtracting the stellar con- 

, tinuum from every galaxy spectrum before making emission line measurements. The 

distribution of the measured Ha Full Width at Half Maxima values of emission line 
galaxies is strongly bimodal, with two populations separated at about 1,200km s~^. 
' This feature provides a natural separation between narrow-line and broad-line AGN. 

The narrow-line AGN are identified using standard emission line ratio diagnostic dia- 
grams. 1,317 broad-line and 3,074 narrow-line AGN are identified from about 100,000 
galaxy spectra selected over 1151 square degrees. This sample is used in a companion 
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, paper to determine the emission-line luminosity function of AGN. 

Subject headings: galaxies: active — galaxies: Seyfert — galaxies: starburst — galaxies: 
^ . quasars: emission lines — surveys 

^Princeton University Observatory, Princeton, NJ 08544 

^Current address; Astronomy Department, Cornell University, Ithaca, NY 14853; liaol@isc.astro.cornell.edu 
^Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721 

^Department of Physics and Astronomy, Johns Hopkins University, 3400 North Charles Street, Baltimore, MD 
21218 

^Max-Planck Institut fiir Astrophysik, D-85748 Garching, Germany 

^Center for Cosmology and Particle Physics, Department of Physics, New York University, 4 Washington Place, 
New York, NY 10003 

^Department of Astronomy and Astrophysics, Pennsylvania State University, University Park, PA 16802 

Apache Point Observatory, P.O. Box 59, Sunspot, NM 88349-0059. 

^Department of Astronomy and National Center for Supercomputer Applications, University of Illinois, 1002 West 
Green Street, Urbana, IL 61801. 



- 2 - 



1. Introduction 

Ever since the definition of Seyfert galaxies (Seyfert 1943) and the first recognition of quasars 
(Schmidt 1963), astronomers have put enormous effort into compiUng large samples of active galac- 
tic nuclei (AGN, in this paper, "AGN" refers to active galactic nuclei at all luminosities, including 

quasars) and trying to understand the physics that powers them. This is not easy, since luminous 
AGN comprise only a few percent of normal galaxies. Based on the distinctive characteristics of 
AGN, different methods have been developed to search for AGN in various wavebands. 

In the optical, AGN show different colors from stars and normal galaxies, especially at high 
luminosity. Many surveys have used color selection (e.g. Schmidt & Green 1983; Boyle et al. 
1990). In particular, the Sloan Digital Sky Survey (SDSS) (York ct al. 2000) uses optical colors to 
identify quasar candidates, which are then observed spectroscopically (Richards et al. 2002). The 
color selection is very efficient but it requires that the optical luminosity of an AGN be at least 
comparable to the luminosity of its host galaxy for the color to be distinctive, and thus the color 
selection systematically misses less luminous AGN at low redshift. 

AGN also show strong optical and ultraviolet emission lines. Broadly speaking, AGN can 
be classified into two types: broad-line and narrow-line AGN. The former show broad permitted 
emission lines, with Full Width at Half Maxima (FWHM) of several thousand km s^"*^, while in 
narrow-line AGN, both the permitted and forbidden emission lines are narrow, with FWHMs ~ 500 
km s~^. This is comparable with emission lines in normal star-forming galaxies, but emission lines 
in narrow-line AGN have considerably greater ionization range. In particular, both high-ionization 
lines such as [NeIII]A3869, [Ne V]A3426 and [OIII]A5007, and low-ionization lines such as [OI]A6300 
and [NI]A5200, are stronger in narrow-line AGN than in normal starforming galaxies. Based on 
this feature, narrow-line AGN can be identified by their distinctive emission line ratios. 

The first line ratio diagram was introduced by Baldwin, Phillips &: Tcrlcvich (1982), who 
suggested that AGN generically have greater [OIII]A5007/H/3 (the flux ratio of [OIII]A5007 to H/5) 
than do galaxies whose emission lines are due to stellar processes. Veilleux & Osterbrock (1987) 
developed their idea and used diagnostic diagrams that consist of combinations of four line ratios: 
[OIII]A5007/H/3, [NII]A6584/Ha, [OI]A6300/Ha and [SII]AA6716, 31/Ha, and further developed a 
semi-empirical line on the diagrams to separate AGN and starburst galaxies. The emission-line 
pairs are chosen specifically so that the two lines in a ratio are at nearly identical wavelengths, 
therefore reddening and spectrophotometric uncertainties are not big effects. These diagnostic 
diagrams have been used ever since as a standard to identify narrow-line AGN. 

Kcwlcy ct al. (2001) developed a set of theoretical separation lines for AGN and starforming 
galaxies on the diagnostic diagrams. By constructing a detailed continuous starburst model with 
large realistic metallicity and ionization parameter ranges, they found that the model folds on the 
diagnostic diagrams and there exist upper limits for starforming galaxies. These upper limits can 
be used to separate AGN and starforming galaxies and they can be fitted to a simple rectangular 



-3- 



hyperbolic shape: 



log 



/[OIII]A5007\ 



f [QIII]A5007^ 



0.61 



+ 1.19 



log 



log([NII]/Ha) - 0.47 
0.72 



+ 1.30 



(1) 



log 



V H/3 ) 
/[OIII]A5007\ 



y R(3 J 



log([SII]/Ha) - 0.32 
0.73 



+ 1.33 



log([OI]/Ha) + 0.59 



These theoretical separation lines between AGN and starforming galaxies are widely accepted 
to use to identify narrow-line AGN in the diagnostic diagrams. 

KaufFmann et al. (2003a), on the other hand, when studying host galaxy properties of narrow- 
line AGN in the SDSS, proposed an empirical and more lenient cut to identify AGN: 



This criterion will select many more galaxies as AGN than Kewley's criteria. In this paper, we 
apply both criteria and discuss their differences in detail. 

There are several spectroscopic galax;y surveys from which people have tried to select AGN 
based on their emission lines. One is the CfA redshift survey (Davis, Huchra & Latham 1983; 
Huchra et al. 1983). Spectra of about 2,400 galaxies were taken to study their large scale distri- 
bution. Huchra, Wyatt & Davis (1982) used these spectra to select AGN by identifying emission 
lines indicative of nonstellar activity. As a result they found roughly 50 Seyfert galaxies, divided 
approximately equally between Seyfert 1 and Seyfert 2 galaxies, making the local Seyfert fraction 
~ 2% (Huchra & Burg 1992). 

Another example is the Revised Shapley-Ames catalog of bright galaxies (Sandage & Tammann 
1981). Among its ~1,300 galaxies, about 50 are Seyfert galaxies of luminosity comparable to those 
found in the CfA Survey. Ho et al. (1997a) took uniform high signal to noise ratio (S/N) spectra of 
486 galaxies using very small apertures centered on the nuclei. 418 galaxies were found to contain 
emission-line nuclei, of which 206 are star- forming, and 211 show AGN activity. This demonstrates 
that low-luminosity AGN activity in the local universe is extremely common, as first discussed by 
Phillips, Charles, & Baldwin (1983). 

Hall et al. (2000) systematically selected high redshift AGN based either on their broad emis- 
sion lines or narrow [NeV] emission lines from the Canadian Network for Observational Cosmology 
field galaxy redshift survey (CN0C2; Yee et al. 2000). They found 47 confirmed and 14 candidate 
AGN in the redshift range 0.27 <z< 4.67. 

The SDSS has opened a new window to the study of AGN by providing a huge number of high 
quality spectra of galaxies. There have been several studies identifying AGN spectroscopically from 
different subsamples of SDSS galaxies (KaufFmann et al. 2003a; Miller et al. 2003). Both studies 
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have focused on naxrow-line AGN identified from the diagnostic diagrams. Miller et al. (2003) 
found AGN signatures in at least 20% of 4,921 galaxies in the redshift range 0.05 < z < 0.095, 
and studied the environment of these AGN. Kauffmann et al. (2003a) selected 22,623 narrow-line 
AGN from a parent sample of 122,808 galaxies and studied their host galaxy properties. The AGN 
detection rate, however, depends very much on the AGN selection criteria used and many other 
details in defining the sample. In this paper, we apply a systematic search for both broad and 
narrow-line AGN in a well-defined sky area in the redshift range < z < 0.33. The framework 
of our narrow- line AGN selection is similar to Kauffmann et al. (2003a), but we differ in many 
details. 

In section 2, we give a brief overview of the SDSS, focusing on those aspects most relevant 
to our study. §3 introduces the parent sample from which we will select our AGN. In §4, we will 
discuss the subtraction of stellar absorption lines from the spectra. The emission line measurements 
are discussed in §5. In §6, the selected AGN are presented and discussed. §7 gives a clean AGN 
sample and we summarize in §8. 

2. The Sloan Digital Sky Survey 

The Sloan Digital Sky Survey (York et al. 2000) is an imaging and spectroscopic survey that 
will eventually cover approximately one-quarter of the Celestial Sphere and collect spectra of ~ 10^ 
galaxies and 10^ quasars. It uses a dedicated 2.5m telescope at Apache Point, New Mexico, with 
a 3 degree field, and a mosaic CCD camera and two fiber-fed double spectrographs to carry out 
the imaging and spectroscopic surveys respectively. A separate 20" photometric telescope is used 
for photometric calibration (Smith et al. 2002; Hogg et al. 2001). The imaging camera (Gunn et 
al. 1998) consists of a mosaic of 30 imaging CCDs with 24 |Um pixels subtending 0.396" on the sky. 
The sky is observed through five broad-band filters {u, g, r, i, z) (Fukugita et al. 1996; Stoughton et 
al. 2002) covering the entire optical band from the atmospheric cutoff in the blue to the sensitivity 
limit of silicon CCDs in the red. The imaging is done in drift-scan mode and the total integration 
time per filter is 54.1s. The 50% completeness limits for point sources are 22.5, 23.2, 22.6, 21.9 
and 20.8 magnitudes respectively and the photometric calibration is reproducible to 3%, 2%, 2%, 
2% and 3% for the five bandpasses, respectively. The image data are processed by a series of 
automated pipelines (Lupton et al. 2001; Stoughton et al. 2002; Pier et al. 2003), which make 
various measurements of the flux of each detected object. 

The SDSS is obtaining spectra of complete samples of three categories of objects: Galaxies, 
Luminous Red Galaxies and Quasars. These spectroscopic targets are selected from the imaging 
data via various target selection criteria. Galaxy target selection is discussed in Strauss et al. (2002). 
Briefly speaking, galaxies are separated from stars by morphology. The magnitude limit cut for the 
galaxy sample was changed several times during commissioning, but is currently r = 17.77, where 
r represents the r band Petrosian magnitude. All magnitudes are corrected for extinction following 
Schlegel et al. (1998). These objects are the sample from which we will select AGN. 
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Quasar target selection (Richards et al. 2002) is based on the nonstellar colors of quasars 
and matching unresolved sources to the FIRST radio catalog. Luminous Red Galaxies are selected 
(Eiscnstein et al. 2001) by a variant of the photometric rcdshift method, aiming to have a uniform, 
approximately volume-limited sample of highly luminous objects with the reddest colors in the rest 
frame to z = 0.5. In this paper, we will not discuss AGN selected from these two samples, unless 
they also satisfy the galaxy sample magnitude cut r = 17.77. 

The SDSS spectra are taken with two fiber-fed spectrographs, covering the wavelength range 
3800-9200 A over 4098 pixels. Each plate can hold 640 fibers, with a fixed aperture of 3". The 
plates are positioned by a tiling algorithm (Blanton et al. 2003) and fibers are assigned to targets. 
Galaxies are among the tiled targets that have the highest priority of having their spectra taken. The 
finite diameter of the fiber cladding prevents fibers on any given plate from being placed closer than 
55" apart. The resolution A/AA varies between 1850 and 2200. The relative spectrophotometry 
is accurate to about 20%. Each spectrum is accompanied by an estimated error per pixel, based 
on photon statistics and the amplitude of sky residuals. The typical S/N for galaxy spectra at the 
sample limit is 16/pixel. 

The spectroscopic data are reduced through the spectroscopic pipelines, spectro2d and specBS. 
Spectro2d reduces the 2-dimensional spectrograms produced by the spectrographs to flux- and 
wavelength-calibrated spectra. SpecBS is different from the SDSS official pipeline, it determines 
classifications and redshifts via a fit to the spectrum in question with a series of rest-frame star, 
galaxy and quasar templates. The basic technique is described by Glazebrook et al. (1998) and 
Bromley et al. (1998). 

Even though specBS has made measurements on emission lines by fitting a single Gaussian at 
positions of each expected emission line, we carry out our own fits in our study. The main reason 
is that a single Gaussian fit is not adequate to model some AGN having both broad and narrow 
emission lines. The emission lines are measured directly from calibrated spectra, as we describe in 
detail in §5. 

3. Parent Sample 

In order to do statistical studies of AGN, it makes sense to start with a complete sample of 
objects within a certain well defined area from the SDSS. We start from 129,625 target objects 
complete in 1151 square degrees. This is about 1/2 of the spectra available in the SDSS Second 
Data Release (Abazajian et al. 2004). Among these objects, there are 98,684 galaxies targeted 
by the main galaxy target algorithm (Strauss et al. 2002; §2) and 17,972 quasar target objects 
(Richards et al. 2002; §2). In this paper, we will focus on the galaxy target objects. However, there 
are 2057 extragalactic objects selected as quasar targets with r(Petrosian) < 17.77, which were 
not targeted as galaxies because they are unresolved, and we include them in our galaxy sample. 
Among these 100,741 galaxies, most of them (99,990) have redshift z < 0.33, guaranteeing that the 
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Ha emission line lies in our spectral coverage. Ha is a very useful emission line in identifying AGN 
(Equation 2), and we limit our galaxy sample to these 99,990 objects. 

4. Stellar Subtraction 

AGN are identified by their emission line characteristics. As described in §2, the SDSS galaxy 
spectra are taken through a fixed 3" aperture, which is large enough to let through not only the light 
from the nucleus but also substantial amounts of stellar light from the host galaxy. For example, 
at the median redshift of the sample (z = 0.1), a 3" aperture subtends about 4h~^ kpc. Moreover, 
galaxies with higher redshift will have a larger host galaxy component in the observed spectra. 
Thus the nuclear emission lines arc often contaminated by the stellar absorption lines of the host 
galaxy. For weak AGN, this contamination can be so severe that the interesting emission lines are 
completely submerged in the absorption lines. Thus before considering AGN selection, we have to 
develop a technique to properly remove the stellar absorption lines. 

The basic idea of stellar subtraction is to build a library of stellar absorption line spectra 
templates, and use them as building blocks to simulate the stellar spectrum of the object in question. 
The library needs to be complete in the sense that it contains enough information on various 
absorption features to be able to simulate the stellar components of various galaxies with widely 
spread metallicities, ages and velocity dispersions. The library is typically composed of star spectra 
generated from a population synthesis code (Bica 1988; Saraiva et al. 2001, Kauffmann ct al. 2003b) 
or direct observational spectra of absorption line galaxies (Ho et al. 1997a) or stars (Engclbracht et 
al. 1998). In this paper the library is constructed by applying the Principal Component Analysis 
(PGA) technique (e.g Connolly et al. 1995; Lahav et al. 1996; Bromley et al. 1998; Eisenstein et 
al. 2003; Yip et al. 2004) to a sample of pure absorption-line galaxies. 

The advantages of building stellar absorption templates via this method are two-fold: only the 
first few eigenspectra are significant, so we can limit the size of the library without losing much 
useful information. Moreover, the eigenspectra are orthogonal to each other, resulting in a unique 
solution to the stellar subtraction fit using these templates. 

4.1. Preparing the Absorption Line Galaxy Sample 

We wish to identify a sample of pure absorption-line galaxies. In practice, we require that the 
Ha equivalent width EW(Ha)<0 (positive EW corresponds to emission lines), and that [OII]A3727 
not be detected. [OII]A3727 is used because it is always apparent even in very weak emission line 
galaxies. Due to the presence of complex absorption near [OII]A3727, we measure this line after 
subtracting a preliminary PCA sample solely with the requirement that EW(Ha)<0. A Gaussian 
function is fit to the residuals between A = 3700 and A = 3754. If the of the fit is less than 
the of a linear fit minus 3, the line is considered significant, and the object is rejected as an 
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absorption-line galaxy. 

By limiting ourselves further to high S/N spectra, we defined three samples of several hundred 
pure absorption line galaxies grouped by their redshifts: 325 galaxies with 0.02 < z < 0.06; 338 
with 0.06 < z < 0.12 and 372 with 0.12 < z < 0.22. PCA is done on each group separately, 
giving three sets of eigenspectra that each can be used to subtract stellar components for galaxies 
of similar redshifts. We divide the sample in three groups in order to obtain a larger wavelength 
coverage for each group as well as the resultant eigenspectra. Figure 1 shows the mean of the 
spectra (i.e., the first eigenspectrum) in each galaxy group. They are essentially identical, although 
the lowest-redshift sample has slightly higher S/N. 

The galaxies within the three groups are shifted to fixed rest-frame wavelength bins using sine 
interpolation. Afterwards, each spectrum is normalized to a constant flux value. Unlike some PCA 
analyses (e.g., Eisenstein et al. 2003), the continua are not subtracted from the spectra before PCA. 
Each sample includes several hundred normalized galaxy spectra, which we express as a matrix S 
of dimensions A'^ x M, where N is the total number of galaxies in each group and M is the total 
number of common wavelength bins. Singular Value Decomposition (SVD) is used to build the 
eigenspectra (Connolly et al. 1995; Bromley et al. 1998). 

4.2. PCA Result and Stellar Subtraction 

The PCA analysis generates a set of eigenspectra, with main features of the absorption-line 
galaxies concentrated in the first few. In this study we will use the first eight eigenspectra as the 
stellar absorption-line templates. Empirical tests show that including more eigenspectra, which 
are basically just noise, does not improve the subtraction further. A minimizing algorithm was 
adopted to determine the synthetic stellar absorption spectrum for each galaxy. The minimizing 
is done over the entire observed wavelength range except regions around the strongest emission 
lines found in AGN and starforming galaxies: Ha, H/?, [NII]AA6548, 84, [OIII]A5007,A4959 and 
[OII]A3727. 

Since the stellar templates are eigenspectra of a sample of pure absorption line galaxies and 
galaxies having young stellar populations tend to have emission lines, the resultant eigenspectra 
mainly represent old stellar spectral features. Thus these eigenspectra are not representative of 
galaxies containing a young stellar component. One possible resolution would be to make sure 
that the PCA sample included enough E+A galaxies, which contain young stellar populations and 
do not have emission lines (Goto et al. 2003; Quintero et al. 2004). In this work, however, we 
simply add an A star spectrum selected from the SDSS spectroscopic data to the absorption-line 
template library to represent the young stellar population. Figure 2 shows an example of the 
stellar subtraction with and without this template, for a galaxy with a young star population. The 
improvement using the A star template is dramatic. 

The stellar subtraction is automatically done to all galaxies in our parent galaxy sample. 
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including those quasars that satisfy the galaxy target selection magnitude cut (§ 2 and 3). However, 
doing stellar subtraction to a bright quasar dominated by a non-thermal continuum in the spectrum 
will certainly be a disaster. We thus add a pure power-law spectrum, written as Cpowcrlaw A^^'^, 
to our template library. Figure 3 shows an example of stellar subtraction for a quasar with and 
without the power-law template. The power-law template significantly improves the subtraction. 
Since the spectrophotometry is not perfect, quasars have a range of intrinsic power-law slopes 
(Richards ct al. 2002), and the power law fit can be systematically in error due to Fell emission 
(Figure 3 and Vanden Berk et al. 2001), the power-law template will certainly not be sufficient 
for clean contimium subtraction for all quasars. However we have found it to be adequate for our 
purposes. For galaxies that do not have a nonthermal power-law component, this template can help 
compensate for continuum shape errors due to errors in spectrophotometry and internal reddening. 

In summary, our template library includes 8 PCA eigenspectra of pure absorption line galaxies, 
a power-law continuum and an A star spectrum. As demonstrated in Figures 2, 3 and further in 
Figure 4, weak emission lines such as H/3 that were originally submerged in the stellar absorption 
line are successfully recovered. Moreover, emission lines such as [OI]A6363 and [NI]A5200 which are 

not apparent at all in the original spectrum clearly stand out after the subtraction. The subtraction 
also helps correct the strength of Ha, which is very important in subsequent AGN identification 
using emission line ratios. 

The stellar subtraction is robust for galaxies of a range of velocity dispersions. Figure 4 
shows two galaxies with very different velocity dispersions (Bernardi et al. 2003); note that in 
both cases, the subtracted spectrum shows no appreciable residuals of the strong absorption lines. 
The eigenspectra include terms that can give absorption lines of different width. This is another 
advantage of the PCA technique over the use of stellar libraries: we do not need to convolve each 
template with a Gaussian broadening function for each galaxy. 

The whole procedure of doing PCA and stellar subtraction using the resultant eigenspectra 
is quite straightforward. Once the stellar component is subtracted, the emission lines can be 
measured, as we now describe. 

5. Emission Line Measurements 

Since we are mainly interested in emission-line galaxies, we will first set up a criterion to remove 
the pure absorption-line galaxies from the parent sample. We require that the EW of the Ha line 
(in the rest frame) be greater than 3 A. In the rare cases in which the Ha line is saturated or 
affected by bad pixels, we examine the equivalent width of the [OIII]A5007 and H/? lines, requiring 
that one of them be greater than 3 A. These equivalent widths are based on the line strength after 
stellar continuum subtraction, divided by the stellar continuum itself. 

This criterion also rejects weak emission line galaxies. This is acceptable since the S/N ratio 
for these lines will be low and thus our ability to distinguish AGN from starbursts will be poor. A 



-9- 



total of 42,435 galaxies, about half of the parent galaxy sample, pass the equivalent width cut. We 
refer to this sample as the emission line galaxy sample. All AGN are selected from this sample. 

The intensity, full width at half maximum (FWHM) , central wavelength and nearby continuum 
value of the main emission lines in these galaxies are measured via Gaussian fits weighted by 
the estimated errors per pixel. The following emission lines are needed for AGN selection: Ha, 
[NII]A6584,48, H/3, [OIII]A5007, [SII]A6716, 31 and [OI]A6300. The following lines are fit with a 
single Gaussian: H/3, [OHI]A5007 and [OI]A6300. The [SH] doublet is fit with a two-Gaussian 
function model. The FWHMs and intensities of the two Gaussians are independent of each other, 
while the central wavelengths are correlated with a single variable and the continuum is shared 
by the two Gaussian functions. Since our fitting is done to a stellar continuum subtracted spectrum 
that has been shifted to rest wavelength, the fitting results for z and the continuum are very close 
to zero. The Ha and [NH]A6548, 84 lines are fitted with three Gaussians. The FWHMs of the [NH] 
lines are kept the same and the intensity ratio of [Nn]A6584 to [Nn]A6548 is fixed to 3, as required 
by the energy level structure of the [NH] ion (Osterbrock 1989). Again, the central wavelengths 
of Ha and the [NH] doublet are correlated with a single redshift parameter and the three lines 
share the same continuum value. The pixels located within 100 A of the central wavelength of the 
emission lines arc used for Gaussian fits (adjacent emission lines in this range are shielded out). 
However, in order to be sensitive to a broad component of Ha, we fit the Ha, [NH] group to the 
range 6565 A ± 300 A (the adjacent [01] and [SH] lines are masked in the fit). We will use the 
of the fit to test for the significance of the broad component, restricting ourselves to the range 
6565 A zb 80 A to reduce the sensitivity to uncertainties in the continuum. 

Some AGN show both broad and narrow permitted emission lines. For these galaxies, a single 
Gaussian function for Ha is obviously not appropriate. To take this into account, we also fit the 

Ha and [NH] doublet with a four-Gaussian model: two for the [NH] lines and two for Ha. The two 
Ha Gaussian functions have the same central wavelength, but difi'erent intensities and FWHMs. 
The final decision of which model (the three-Gaussian model or the four-Gaussian model) to use 
for a given galaxy is done by comparing the of the two model fits, xi and xl respectively. In 
particular, we choose the four-Gaussian model fit when: 

(xi-xl-2)/xl > 0.2 (3) 

This criterion is empirical, but is inspired by a similar statistic for linear fitting models (Lupton 
1993). If we would like to fit a data set {xi,yi) with a linear fitting model ymodel = y{x)^ and the 
errors aj are Gaussian, then 



X 



2 _ I Hi ymodel \ 



2 



follows a Xn-A; distribution, where n and k are the numbers of the data points and the parameters 
used in the model respectively. If there are two linear models M and TV, with k and {k — r) 
parameters respectively (i.e. model M has r extra parameters), their x^ functions Xatj Xm ^^^^ 
then follow Xn-fe X^_(jfc_r) distributions. It can be proved that for linear models, Xm and 
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Xm ~ Xn independent, thus the quantity 

^ " Xh/in - k) 

follows a Fr^n-k distribution. If / is large, then we can accept the hypothesis that the added 
parameters significantly improve the fit. 

In our case, our models are non-linear in the parameters, and Xm ^^'^ Xm ~ Xn will not 
be independent. The four-Gaussian model has two more degrees of freedom than does the three- 
Gaussian model, so the numerator of the criterion is written as (x| — xl — 2). The limit "0.2" 
in equation (3) is an empirical number and is demonstrated to be appropriate from our manual 
inspection. 

However it is not always unambiguous to identify the broad Ha component using the above 
criterion, since not all emission lines are well-fit with Gaussians (Strateva et al. 2003). In particular, 
narrow emission lines generally have extended wings at their bases (Ho et al. 1997b). If the narrow 
emission line is strong, this non-Gaussian feature will become prominent, and a 4-Gaussian model 
will be chosen by Equation 3. In this case, the height of the broader component hi will be small 
compared to the height of the narrow component, ^2 and the Gaussian width of the broader 
component a will be relatively small. To not count such cases as broad line, we stick with the 
3-Gaussian fit for those objects which satisfy, 

cr < 20A(~ 2200km s~^) and /n//i2 < 0.1 (6) 

no matter what the criterion of equation (3) might indicate. 

We compared our empirical criterion (equation (6)) with the well defined Bayesian Information 
Criterion (BIG, Liddle 2004). We found that we are as efficient in choosing broad- line component 
as BIG, and at the same time, our criterion identifies far fewer fake broad-line components than 
BIG docs. 

After fitting all relevant emission lines, the line strengths are calculated based on the line fitting 
parameters. The observed SDSS spectra are accompanied by an estimated error per pixel, based on 
photon and read noise statistics and variation among sky spectra. These errors are approximately 
independent between pixels and are good to 8% (McDonald et al. 2004). In stellar continuum 
subtraction wc assumed the stellar templates perfect since they are very high S/N. The errors in 
the line strengths are calculated using standard propagation of errors. Those emission lines that 
are less than 3a detections are marked as weak. If they are needed in identifying an AGN, special 
care should be applied (§6.2). 
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6. AGN Statistics 

6.1. Broad-Line AGN 

The measured emission line parameters can be used to identify AGN from the emission hne 
galaxy sample. Wc select broad-line AGN by checking the Ha line width. Figure 5 gives the 
distribution of the FWHM values of the Ha emission line for over 40,000 emission line galaxies. 
For galaxies that prefer two Gaussian components for Ha (§5), the width of the broader component 
is plotted. The distribution is clearly bimodal, with a minimum at 1,200 km s~^. The first peak is 
at about 200 km s~^, which reflects the resolution limit of the SDSS spectrographs. It extends far 
above the limit of the plot, demonstrating that most of the galaxies in the sample are starforming 
galaxies with narrow Ha lines. Clearly separated from these normal galaxies, the second group of 
galaxies has substantially larger Ha FWHM value. Naturally, we use FWHM(Ha)> 1,200km 
as the selection criterion for defining broad-line AGN. 

The inserted plot in Figure 5 shows the Ha FWHM distribution for AGN only (the detailed 
selection of narrow-line AGN is described in §6.2). Narrow-line AGN have typical Ha FWHM 
values similar to those of usual emission-line galaxies. The paucity of objects with FWHM ~ 1, 200 
km s"^ is not well understood and will need to be explained in AGN unification models. 

However, we should make sure that the bimodality is not a consequence of our fitting procedure 
in choosing the 4-Gaussian fitting model over the 3-Gaussian model. To do so, we intentionally 
change the limiting numbers of the criteria described in §5. We write the criteria (equations (3) 
and (6)) as: 

(xi - Xl - 2)/xl > A and {a > B or h/h2 > C) (7) 

Our default values arc A = 0.2, B = 20A (2182 km s'^), C = 0.1. Now we change this set of 
numbers to ^ = 0.6, B = 30A (3249 km s^^), C = 0.7 separately. The changes will all reduce 
the number of objects with broad Ha components. However, those galaxies which are distinctively 
broad will still be selected by the new criteria. The top two panels of Figure 6 shows the Ha FWHM 
distribution with the changed criteria. As can be seen, no matter how the criteria are changed, 
the bimodality feature and the minimum at about 1,200 km remain the same even though the 
number of selected broad-line AGN changes substantially. A complementary test is to loosen the 
criteria a little. Out of the concern that the B = 20k cut might artificially reject objects with 
FWHM around 1,200 km s^^, we change this criterion to -B = oA. The bottom panel of Figure 6 
shows the Ha FWHM distribution with A = 0.2, B = A and C = 0, 0.1, 0.7 respectively. This 
yields many false broad-line AGN. However, the bimodality and the location of the minimum are 
still the same. Therefore, the bimodal feature is insensitive to the specific details of the criteria in 
equation (7). 

One might also suspect that the bimodality is caused by the fact that the values we have used 
in the plot are coming from two different fitting models. In the distribution plot, the values of the 
narrower group usually come from the fitting parameters of the 3-Gaussian model while those of 
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the broader group are mostly from the 4-Gaussian model. The change from one model to another 

might cause some discontinuity in an otherwise continuous distribution. To test this, Figure 7 
shows the distribution of width using the 3-Gaussian model fits alone; the bimodal feature is still 
apparent. 

As a final check of the reality of the deficit of galaxies with FWHM~ 1, 200 km s~^, we carried 
out simulations of galaxies with composite Ha profiles, in which the broad component had 1200 
km s~^ < FWHM < 2200 km s~^ and 0.1 < /ii//i2 < 0.3, and to which we add a stellar continuum 
and characteristic noise. After putting these objects through our code, we found that we would 
miss only a small fraction of broad-line AGN with Ha FWHM close to 1,200 km s"^, not nearly 
enough to explain the bimodality. Therefore, the bimodal feature is clearly real. 

Defining broad-line AGN as objects with Ha FWHM > 1,200km s~\ 1317 broad-line AGN 

are identified from 42,435 emission line galaxies. For those galaxies that need two Ha Gaussian 
functions but are not classified as broad (FWHM(Ha)< 1,200km s^-^), wc use a 3-Gaussian fit for 
Ha and [Nil] so that we can have consistent Ha measurements when they are used in emission line 
ratio diagnostic diagrams, as we now describe. 

6.2. Narrow-Line AGN 

Narrow-line AGN cannot be separated cleanly from star-forming galaxies by their emission line 
widths. Therefore, we will use the traditional emission line ratio diagnostic diagrams (Veilleux & 
Osterbrock 1987; Section 1) to select them. In Figure 8 we place the galaxies in the emission line 
galaxy sample in these diagrams (broad- line AGN selected in Section 6.1 are excluded). Galaxies 
with relevant emission lines detected with less than 3a significance are omitted in the diagrams. 
There are roughly 40,000 galaxies contained in the [Nil] / Ha and [SII]/ Ha diagrams, but only about 
23,000 galaxies in the [01]/ Ha diagram. Due to the large number of galaxies available, we can see 
a continuous distribution in the diagrams. AGN have larger [0III]/H/3, [Nil]/ Ha, [01]/ Ha and 
[SII]/ Ha values than do star bursts, and are clearly separated from the loci formed by starforming 
galaxies. The AGN and starforming galaxy separation lines developed by Kaufi'mann et al. (2003a, 
short-dashed), Kewley et al. (2001, solid) and Veilleux & Osterbrock (1987, dash-dot) arc plotted 
in the diagram. Kauffmann's and Kewley's criteria agree much better with the shape of the density 
profile; the disagreement between Kewley's lines and the data are almost all in the iO.ldex error 
range (dashed lines). 

Within Kewley's separation scheme, we can identify AGN as those located to the upper right of 
the lines. About 10% of all emission line galaxies are classified differently in the three diagrams. We 
will refer to these as mismatched galaxies. One half of these are classified as AGN by the [01]/ Ha 

diagram, but not by the other two. The same discrepancy was also noticed by Stasihska & Leitherer 
(1996) and Dopita et al. (2000), who ascribed the enhancement of [01] lines in starburst galaxies 
to the mechanical energy released into the gas by supernovae and stellar winds. Dopita (1997) 
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noted that when even a weak shock compresses the gas, the postshock local ionization is sharply 
lowered, leading to a strong enhancement of the [01] lines. Another quarter of the mismatched 
galaxies arc those classified as AGN in the [SIX]/ Ha diagram but not in the [Nil] /Ha diagram, 
and the [01] line is of too low significance to use the [01]/ Ha diagram. The cause could again be 
shock excitation: the relatively cool high-density regions formed behind the shock front emit strong 
[Sn] emission lines (Dopita et al. 2000). The remaining quarter of the mismatched galaxies are 
randomly distributed, possibly due to noise. Overall, about 95% of the mismatched objects have 
enhanced [01] or [SIX] emission lines compared to the [NXX] lines, perhaps due to shock excitation 
from supernovae and stellar winds in a starburst galaxy. The small overall number of mismatches 
demonstrates that Kewley's theoretical separation lines give very consistent excitation mechanism 
classifications among the three diagrams. 

Using Kewley's scheme, we define the galaxies located to the upper right of the lines in all 
three diagrams as narrow-line AGN. If the emission lines used in the diagrams are weak (§5), we 
take special care: if the [Nil] or [OIII] line is found to be weak, the galaxy will not be identified as 
an AGN unless the Ha line is also broad. However, if the [01] line or [SII] line is found to be weak, 
then the classification based on the [01]/ Ha or [SII]/ Ha diagram is defined to be inconclusive 
and the identification of the galaxy depends on the remaining diagrams. Galaxies with one or two 
inconclusive classifications and AGN classifications in the remaining diagram(s) are also included 
in our sample. We do not include any mismatched objects in the final narrow-line AGN sample. 
Based on these criteria, 3074 narrow-line AGN are selected. 

The selection criteria are different from those used by Kauffmann et al. (2003a) (Equation 2, 
Figure 8 short-dashed line). In the [Nil]/ Ha vs. [0III]/H/3 diagram, galaxies are distributed 
on two arms. Normal galaxies have weaker [Nil] /Ha ratio than do AGN. So it is expected that 
they form the locus on the left arm (with smaller [Nil] / Ha values) and galaxies on the right arm 
are AGN. But Kewley's line cuts through the right arm. Kauffmann et al. propose to classify all 
galaxies located on the right arm as AGN, generating an empirical line separating the two branches. 
This line lies well below Kewley's criterion in the [Nil]/ Ha diagram and the number of galaxies 
which fall between the two criteria is very large. 

This difference between the two criteria will cause tremendous disagreement on AGN statistics. 
If Kauffmann's criterion is adopted, we would have selected about 10,700 narrow-line AGN, roughly 
three times the AGN selected via Kewley's criterion. To understand the disagreement better, we 
divide the galaxies in the right arm into several groups by their locations along the arm; each group 
is assigned a different color in Figure 9. These galaxies are accordingly plotted in the [01] /Ha 
and [SII]/ Ha vs. [0III]/H/3 diagrams. In addition, we plot the Ha FWHM distribution for the 
galaxies in each group. Very interestingly, galaxies in the right arm form a similar sequence in the 
three diagrams: galaxies with larger [Nil]/ Ha values usually have larger [01]/ Ha and [SII]/ Ha 
values. Furthermore, galaxies with larger [Nil]/ Ha or [0III]/H/3 values typically have broader Ha 
emission lines. This gives us a hint that some galaxies on the right arm might be galaxies with 
mixed AGN and starburst components, with a larger AGN component corresponding to a higher 
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position on the diagrams (Kauffmann et al. 2003a). However, some of the AGN and starburst 

composite galaxies are degenerate with starforming galaxies with widely spread metallicity and 
ionization parameters, especially in the [01]/ Ha and [SHJ/Ha vs. [Oni]/H/3 diagrams. Because 
of the degeneracy, Kewley's separation lines, which are the upper limit for pure starburst galaxies, 
are unable to identify every galaxy with an AGN contribution. 

Therefore, the AGN sample created using Kauffmann's criterion includes these AGN + star- 
burst composite galaxies; they make up the majority of the sample. While these galaxies include an 
appreciable contamination from star formation in their emission line strengths, the star formation 
contribution to [OHI] is relatively weak (Kauffmann et al. 2003a) . The AGN sample created using 
Kewley's criterion includes only those that are unambiguously AGN dominated. Both samples have 
their pros and cons, so we will keep both samples. In Hao et al. (2005; paper H), the luminosity 
functions are measured for both samples and the results are compared. 

Some of the narrow-line AGN have very strong low ionization lines; these objects are "Low 
Ionization Nuclear Emission Line Regions" ( "LINERs" ) . Ever since the first definition of the LINER 
class by Heckman (1980), the excitation mechanism causing LINERs has remained unknown. There 
is vigorous debate on whether or not LINERs belong to the AGN family. Filippenko et al. (1985) 
showed that many LINERs have a broad emission-line component, suggesting they are indeed AGN. 
Kewley et al. (2001) established a series of starburst models and developed a set of extreme mixing 
lines, below which the galaxies can only be modeled by pure shocks without an ionizing precursor, or 
a power-law ionizing radiation field with an extremely low ionization parameter. The long-dashed 
lines in Figure 8 show the extreme mixing line in the [01]/ Ha diagram. If this line alone is adopted 
here in selecting LINERs, there are 650 LINERs identified in the narrow-line AGN sample. These 
objects will be studied in detail in future work and will not be discussed further here. 

It is interesting to know how the narrow components of broad-line AGN are located on the 
diagrams. In Figure 10, we select those broad-line AGN with real narrow Ha components (whose 
line widths arc less than 5 A) and plot their emission line ratio with the narrow components of 
Ha and H/5 in the diagram. Most of the broad-line AGN are concentrated on the right arm 
of the [Nil]/ Ha vs. [0III]/H/3 diagram, indicating a strong AGN contribution. Therefore, the 
classifications of these galaxies from the narrow component and the broad component are consistent. 
There are also some broad-line AGN with their narrow components corresponding to locations to 
the right of the mixing line in the [01]/ Ha diagram (Figure 8), these can be characterized as 
LINERs (Filippenko et al. 1985). 

The final AGN sample includes 1,317 broad-line AGN, 3,074 narrow-line AGN if Kewley's 
criteria are adopted and 10,700 narrow-line AGN with Kauff'mann's criteria. This AGN sample is 
complete over 1151 square degrees, defined by the criteria described in this paper^°. 



The sample is available at: littp://isc. astro. cornell.edu/~haol/agn/agncatalogue.txt 
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7. Summary 

The large number of high quahty spectra available in the SDSS provides us with great oppor- 
tunities for AGN studies. In this paper, we have systematically identified AGN from the SDSS 
spectroscopic data from their emission line characteristics. In order to make unbiased measure- 
ments of the emission lines, we have carefully subtracted absorption stellar spectra from the ob- 
served spectra. The absorption line templates are constructed by applying PCA to pure absorption 
line galaxies. Emission lines are measured from corrected emission line spectra by fitting Gaussian 
functions. In particular, to obtain proper identification of broad emission lines, the Ha and [Nil] 
emission lines are fitted with both a three-Gaussian model and a four-Gaussian model. The final 
decision of which model to use is based on a comparison of values of the two model fits. 

Using measured parameters of the Ha emission lines, broad-line AGN are identified. It is 
found that the distribution of Ha FWHM values of emission-line galaxies in the parent sample is 
bimodal, with two populations separated at about 1,200km s~^. Various tests have been applied 
to confirm that this feature is not caused by selection effects. The bimodal feature also gives a 
natural separation between broad-line and narrow-line AGN, and is significant in understanding 
the "unified theory" of AGN. A naive explanation could be that the physical parameters of the 
torus are correlated with the black hole mass. For an AGN with a black hole mass less than a 
certain value, the dust torus becomes so large that it almost totally blocks the broad-line region. 
Therefore, it is hard to observe any broad-line component with small Ha FWHM values. 

Narrow-line AGN are identified using traditional diagnostic diagrams. Due to the large number 
of galaxies available, starforming galaxies follow well-defined loci. Theoretical separation lines 
between AGN and starforming galaxies developed by Kewley et al. (2001) agree fairly well with 
the shape of the galaxy density contour and bring few mismatched classifications among the three 
diagrams. Only in the [Nil]/ Ha diagram does Kewley 's line seem to disagree with the data, where 
the galaxies are distributed into two arms and Kewley's line cuts through the right arm (larger 
[Nil]/ Ha value). Kauffmann et al. (2003a) argued that all objects in the right arm should be 
considered as AGN. Based on the SDSS data, they proposed to use the empirical line separating 
the two arms to select AGN. This line is lower than Kewley's line, and would classify many more 
objects than Kewley's criterion as narrow-line AGN. We demonstrated that AGN selected via 
Kewley's criterion only include those that are dominated by active nuclear activity, while AGN 
selected via Kauffmann's criterion also includes AGN -|- starburst composite galaxies. Both samples 
have their pros and cons, and we keep both for further studies. 

The final AGN sample, complete to magnitude cut r(Petrosian) < 17.77, includes 1317 broad- 
line and 3074 narrow-line AGN (10,700 if Kauffmann's criterion is used) over 1151 square degrees. 
This is based on 1/8 the eventual SDSS survey; the sample will grow substantially. This sample is 

very useful for various studies. First, with the large number of AGN in this sample, we can evaluate 
the AGN luminosity functions; we carry out this analysis in Paper II. Second, by correlating these 
AGN with their host galaxy parameters, we can study the relationship between AGN and their 
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host galaxies for both narrow-line (Kauffmann et al. 2003a) and broad-line AGN. Third, from 

the velocity dispersion of the AGN host galaxies, we can directly measure the approximate black 
hole mass using the Mbh — o" relation developed first in quiescent galaxies (Gebhardt et al. 2000; 
Ferrarese &; Merritt 2000; Ferrarese et al. 2001; Tremaine et al. 2002), and further derive the 
accretion rate for various types of AGN (Heckman et al. 2004). This will help us to understand 
their accretion mechanisms. Fourth, we can measure the cluster properties and local environments 
of different types of AGN in the sample (Miller et al. 2003), which is important to constrain AGN 
formation scenario. Furthermore, we can also correlate this AGN sample with surveys in other 
wavebands, such as X-ray, radio, infrared, etc. to better understand the physics of AGN. Studies 
on these topics will be covered in future papers. 
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Fig. 1. — The mean spectra of pure absorption-line galaxies in three groups of redshift. Within the 
common wavelength range, the mean spectra are almost identical, although the S/N is highest at 
lowest redshift. PCA is applied to each group and three sets of eigenspectra are obtained. 
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Fig. 2. — A subtraction example (a normal starforming galaxy) comparing stellar-subtraction 
without (top) and with (bottom) an A star spectrum included among the absorption-line templates. 
The grey spectra are the stellar spectra constructed using the absorption-line templates. The A 
star spectrum significantly improves the performance of the stellar subtraction in this case. Note 
that the subtraction makes weak emission lines such as [OI]A6363 more prominent. 
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Fig. 3. — A subtraction example comparing stellar-subtraction without (top) and with (bottom) a 
power-law spectrum included in the absorption-line templates. The power-law spectrum helps to 
properly subtract the quasar continuum without generating fake emission lines. Note the presence 
of Fell A4570 here, which does affect the continuum fit. 
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Fig. 4. — The subtraction using the same set of eigenspectra templates works well for galaxies with 
very different host stellar velocity dispersions. The color notation is the same as in Figure 2. Note 
in particular that there are no substantial residuals around the strong absorption lines in the two 
cases. 
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Fig. 5. — The Ha FWHM distribution for emission line galaxies. The inserted plot is the Ha 
FWHM distribution for narrow-line AGN and broad-line AGN after removing star-forming galaxies. 
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Fig. 6. — The distribution of Ha width for emission line galaxies with various criteria (cf, equa- 
tion (7)) for choosing the four-Gaussian model for the Ha and [NH] lines over the three-Gaussian 
model. The existence of the bimodal feature is insensitive to the exact criteria used. 
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Fig. 7. — The Ha FWHM distribution of emission line galaxies, when the Ha and [NH] lines are 
fitted with a three-Gaussian model in all galaxies. 
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Fig. 8. — Emission-line diagnostic diagrams with the separation hnes taken from Kewley et al. 
(2001) (soHd), Kauffmann et al. (2003a) (short-dashed) and Veilleux & Osterbrock (1987) (dot- 
dashed). The dotted lines are the ±0.1 dex of Kewley 's separation. The contour plots are density 
contours of narrow emission line galaxies with relevant emission lines detected at least 3(t signif- 
icance (the typical error is shown in the lower left corner). The long-dashed line in the [01]/ Ha 
vs. [0III]/H/3 diagram is the extreme mixing line (Kewley et al. 2001), to the right of which are 
possible AGN with very low ionizations, i.e. LINERs. 
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Fig. 9. — Diagnostic diagrams with Kewley's separation lines (solid). The galaxies are colored by 
their locations on the [Nil] /Ha diagram. In the lower right figure, the distribution of Ha FWHM 
of galaxies in each color group is plotted in the corresponding color. The typical error is shown in 
the lower left corner. 
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Fig. 10. — Diagnostic diagrams with large points representing the locations of narrow components 
of broad-line AGN. The solid lines, short dashed line and errorbars are the same as in Figure 8. 
Note that most points are located in the AGN region of the diagram. 



